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The primordial perturbation is widely accepted to be generated through the vacuum fluctuation of the scalar 
field which drives inflation. It is, however, not completely clear what is the natural vacuum in the inflationary 
universe particularly in the presence of non-linear interactions. In this series of papers, we will address this 
issue, focusing on the condition required for the removal of the divergence from the infrared contribution to 
loop diagrams. The IR regularity is deeply related to preserving the coordinate independence in describing the 
- - - local observable universe. In our previous works, the IR regularity condition was discussed using the slow roll 

I expansion, which restricts the background evolution of the inflationary universe. We will show more generally 

. that the absence of infrared divergence leads to non-trivial constraints on the allowed quantum states. 

o; 

C/^ ■ I. INTRODUCTION 

o" 

T— I ■ It is widely accepted that primordial curvature perturbations originate from the vacuum fluctuation of the inflaton field in the 

■ inflationary universe. Different choices of quantum states lead to different statistics of the primordial fluctuation. Among various 

(— I ■ options of vacuum states, we usually choose the adiabatic vacuum. A free scalar field can be understood as a set of independent 
» , harmonic oscillators. Taking the adiabatic vacuum looks reasonable, because each oscillator tends to have a fixed frequency as 

■ the time scale becomes shorter than that of the cosmic expansion. However, once we start to take into account the self-interaction 
1^ of the field, the reasoning of choosing the adiabatic vacuum becomes obscure. Therefore it is an interesting question whether 

• some physical requirement prohibits taking an arbitrary quantum state or not. This is the question we address in this paper 
] The key idea is to focus on fluctuations that we can measure in observations. Since the region we can observe is limited to a 
. portion of the whole universe, we need to take into account the local property of the observable fluctuations. In case we can get 
J> ■ access to information only on a portion of the whole universe, there appear additional degrees of freedom to choose boundary 
■"sj" . conditions in fixing the coordinates. As shown in our recent publications [1-6] and will be briefly described later, these degrees 

■ of freedom in the boundary conditions can be thought of as residual degrees of freedom in fixing the coordinates of the local 
] observable universe. The observable fluctuations should be insensitive to the coordinate choice in the local observable universe. 

I • We showed that removing this ambiguity significantly affects the infrared (IR) behaviour of primordial fluctuations. In the 
On ] conventional gauge-invariant perturbation, in which this ambiguity is not taken care of, it is widely known that loop corrections 
. of massless fields such as the inflaton diverge because of the IR contributions ll7l- l28ll . One may expect that performing the 
quantization after we remove these residual coordinate degrees of freedom in the local universe can cure the singular behaviour 
of the IR contributions. However, these residual coordinate degrees of freedom are not present if we deal with the whole universe 
including the spatial infinity, since what we call the residual coordinate transformation diverges at spatial infinity. In this sense, 
if we think of the whole universe, they are not gauge degrees of freedom which preserve the action to be invariant. They 
can be removed only when our concern is limited to our local observable region of the universe. However, if we consider to 
[ quantize the system excluding the region outside of our observable universe from the beginning, the selected quantum state, in 
general, depends on the choice of the boundary conditions. In Ref. [[Ij], we assumed that the boundary conditions are imposed 
at a finite distance to remove the degrees of freedom in choosing coordinates, but the quantization (or equivalently setting the 
initial quantum sate) was performed considering the whole spatial section of the universe. In that analysis we concluded that 
there is no IR divergences irrespective of the initial quantum state. However, there were several overlooked points. After we 
set the initial quantum state at a finite initial time, we need to translate it into the language written solely in terms of the local 
quantities. Namely, the Heisenberg operators corresponding to the perturbation variables should be transformed into the ones 
that satisfy given boundary conditions using an appropriate residual spatial coordinate transformation. Since the non-linear part 
of this transformation is not guaranteed to be IR regular, it can produce IR divergences. 

If we send the initial time to the past infinity, one may think that this problem may disappear because our observable region 
in comoving coordinates becomes infinitely large in this limit. However, this argument is not so obvious. If we think of, for 
instance, the conformal diagram of de Sitter space, the causal past of an observer never covers the whole region of a time constant 
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slice in flat chart. Furthermore, in this case we also need to care about the secular growth of fluctuation. In Ref. iQl], we discussed 
the absence of the secular growth, focusing only on modes beyond the horizon length scale. To complete the discussion, we 
also need to include the vertexes which contain modes far below the horizon length scale, because these vertexes may yield the 
secular growth through correlations between sub horizon modes and the super horizon modes. 

Later in Refs. leaving the residual coordinate degrees of freedom unfixed, we calculated quantities that are invariant 

under the residual coordinate transformation in the local universe under the so called slow-roll approximation. Then, we derived 
the conditions on the initial states that guarantee the IR regularity. These conditions are the ones that guarantee the absence of 
additional IR divergences originating from the coordinate transformation mentioned above. However, the physical meaning of 
these conditions is not transparent from their derivation in Refs. iHQ]. 

In this paper, as in the case of the conventional perturbation theory and also as in Refs. IHQl, we will perform the quantization 
over the whole universe. Then, calculating observable quantities, which are invariant under the residual coordinate transforma- 
tion in the local universe, we will derive the necessary conditions for the initial quantum state to be free from IR divergences 
more generally, removing the limitation to the slow-roll approximation. We will also show that the IR regularity conditions can 
be thought of as the conditions that request the equivalence between two systems related by means of a sort of scale transforma- 
tion of spatial coordinates. In this paper, we will consider a case in which the interaction is turned on at a finite initial time. In 
the next paper, we will discuss the case when the initial time is sent to the past infinity. This paper is organized as follows. In 
Sec. [n] we first show a symmetry of the system, which plays an important role in our discussion. Then, following Ref. lUlHl, we 
construct an operator which remains invariant under the residual coordinate transformation. In Sec.|IIIl we address the conditions 
on the initial states. In Sec.|IV] we will summarize the results of this paper and will outline future issues. 



II. PREPARATION 



In this section, we will introduce basic ingredients to calculate observable fluctuations. In this paper, we consider a standard 
single field inflation model whose action takes the form 

S=^JV^[R- - 2Vi<t>)]d'^ , (2.1) 

where Mpi is the Planck mass and (f> is a dimensional scalar field divided by Mp\. In Sec. Ill Al we show a symmetry of this 
system, which plays a key role in our argument. In Sec. IIIBI after we introduce a variable which preserves invariance regarding 
the coordinate choice in the local universe, we provide a way of quantizing the system. 



A. Symmetry of the system 



To fix the time slicing, we adopt the uniform field gauge Scf) ~ 0. Under the ADM metric decomposition, which is given by 

ds^ = -N^dt^ + h,j{dx' + N'dt){dx^ + N^dt) , (2.2) 

we take the spatial metric hij as 

h,^=e'iP+C)[,s,^^^ , (2.3) 
where denotes the background scale factor, ( is the so-called curvature perturbation and 5^ij is traceless: 

S-f', = , (2.4) 
where the spatial index was raised by (5'^ . As spatial gauge conditions we impose the transverse conditions on S%j : 

d,Sf, = . (2.5) 

In this paper, we neglect the vector and tensor perturbations. The tensor perturbation, which is a massless field, can also 
contribute to the IR divergence of loop corrections. We will address this issue in our future publication. 

Solving the Hamiltonian and momentum constraints, we can derive the action which is expressed only in terms of the curvature 
perturbation 1I29I1 . Since the spatial metric is given in the form e^'^e^'^dx^, we naively expect that the explicit dependence on the 
spatial coordinates appear in the action for ( only in the form of the physical distance ePe'^dx. We first examine this property. 

Using the Lagrangian density in the physical coordinates Lphys, we express the action as 



5= / dtdxL{x)^ / dtdxe^^P+'^^Lphysix) 



(2.6) 
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We can confirm that the Lagrangian density i^phys is composed of terms which are in a covariant form regarding a spatial 
coordinate transformation such as h^^diNj and h^^ Nidj(. Using 

dx' = eP+'^dx' , d/dx' = e-'^P+'^^d/dx' , N, = e^^^+^^iVi , (2.7) 
we can absorb the curvature perturbation ( without differentiation which appears from the spatial metric, for instance, as 

h^^m, = S'^^e-'^P+Od^N^ ^ {J_N^ + N,^^^ . (2.8) 
After a straightforward repetition, the Lagrangian density can be recast into 

5 = J dtdxe^'-P+i^CphyslDCN, N,]. (2.9) 

Here, to stress the fact that all which remain unabsorbed are associated with differentiation dt or d/dx^, we introduced T) 
which denotes differentiations in general. Using this expression of the action, the Hamiltonian constraint is given by 

e[VCN, m := = e-^^^^^)^ = 0, (2.10) 



and the momentum constraints are given by 



e.pCiV, m := ^^e^ = e-2(p+';)^ = 0. (2.11) 

Perturbing these constraints, we can show that the Lagrange multipliers N and Ni are given by solutions of the elliptic type 
equations: 

^) iv = /[a)C(a^)], N, = !m{^)\. (2.12) 

To address the system of the whole universe, the integration in the action is assumed to be taken over the whole universe. Then, 
by assuming the regularity at the spatial infinity, these Poisson equations can be uniquely solved as 

iV = iV[2)C(x)] , TV, = 7Y,[DC(x)]. (2.13) 

Substituting these expressions of and Ni into the action ( 12.91 ), we can show that the action takes the form: 

dtd^xliphys{Ti^^C{x)\. (2.14) 

Here, we put the subscript da; on T) to specify the spatial coordinates used in CD. The equation ( |2.14| i explicitly shows that the 
action for only depends on the physical spatial distance. 

Now we are ready to show the dilatation symmetry of the system, which plays a crucial role in discussing the IR regularity. 
(See also the discussions in i30-36ij.) Changing the coordinates in the integral from x to e^^x with a constant parameter s, we 
can rewrite the action (12.14b as 

^ dtd^xe^^P^'-^LphA^^.^.^^Cixy^ = y dtd3a;e3(''+^-^)£p„,,[2),p+c-.dxC(i, e'^a;)] . (2.15) 

This implies that the action for C possesses the dilatation symmetry 

5" = ^ dtd?xli{C,[x)\ ^ j dtd^xL [Q{t, er'x)-s\ . (2.16) 

As long as we consider a theory which preserves the three-dimensional diffeomorphism invariance, the system preserves the 
dilatation symmetry. It is because the above-mentioned dilatation symmetry can be thought of as one of the spatial coordinate 
transformations. 

To make use of the dilatation symmetry, we introduce another set of canonical variables than (^{x) and its conjugate momentum 
tt{x). We can show that 

Cix) := Cit, e-'x) , ji{x) := e"^"7r(t, e^^x) (2.17) 
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satisfy the canonical commutation relations as well as ({x) and tt{x). Actually, using the commutation relations for ({x) and 
7r(.T), we can verify 



and also 



at, x), n{t, y)] = [C(t, e-'x), er^'^it, e'^y)] = e-^^zJ^^) (^-^(a; - y)) = z5(3)(a, - y) , (2.18) 



C(i, a;), C(t, y) ^[^{t,x),^{t,y)]=Q. (2.19) 



Using Eq. ( |2.16t . we can show that the Hamiltonian densities expressed in terms of these two sets of the canonical variables 
are related with each other as 

d?x'K[m,<x)] = f d'x[n{x)ax)-m^)]} 

d^x |e-3''7r(t, e-^xXit, e'-'x) - £[C(t, e'^'x) - s]| 

= J d^x^[C{x) - s, n{x)] =: J d^xJ<[C{x), n{x)] , (2.20) 

where on the second equahty, we again changed the coordinates in the integral as a; e^^a;. Note that the Hamiltonian density 
for and tt is given by the same functional as the one for ( and tt with ( shifted hy ( — s. 

We define the non-interacting part of the Hamiltonian density for ( and tt by the quadratic part of J{ in perturbation, assuming 
that s is as small as C, and tt , as follows. 



*Co[C», nix)]^5{„[ax)~s, iT{x)] 



(2.21) 



where J£o is the free part of the Hamiltonian density for ( and tt. Using Eq. (12.20b with Eq. ( I2.21l i. we can show that the 
interaction Hamiltonian densities which are defined by Jij := J{ — J{o for each set of the canonical variables are related as 

5f/[C(x), nix)] = :Ki[C{x) - s, nix)] . (2.22) 

Note that the Hamiltonian densities Jf / and !K/ also take the same functional form except for the constant shift of (. 



B. Residual coordinate transformations and quantization 

In this subsection, we consider a way to calculate observable fluctuations. First, we begin with the classical theory. To obtain 
the action for ( by eliminating the Lagrange multipliers N and iVj, we need to solve the Hamiltonian and momentum constraint 
equations. As is schematically expressed in Eqs. (12.12b . these constraint equations are given by the elliptic-type equations. Here 
we note that the region where we can observationally access is restricted to the causally connected region whose spatial volume 
is bounded at a finite past. As far as we are concerned only with the observable region, boundary conditions for Eqs. (12.12b 
cannot be restricted from the regularity of the spatial infinity, which is far outside of the observable region. The degrees of 
freedom for solutions of N and Ni can be understood as degrees of freedom in choosing coordinates. Since the time slicing 
is fixed by the gauge condition (50 = 0, there are remaining degrees of freedom only in choosing the spatial coordinates. In 
the following, we refer to such degrees of freedom which cannot be uniquely specified without the knowledge from the outside 
of the observable region as the residual "gauge" degrees of freedom in the local universe. We put " ", because changing the 
boundary conditions for N and Ni in the local region modifies the action for C, obtained after solving constraint equations, while 
it keeps the action expressed by N, Ni, and ( invariant. In this sense, the change of the boundary condition is distinct from the 
usual gauge transformation, which keeps the action invariant. 

The observable fluctuations should be free from such residual gauge degrees of freedom. Following Refs. IHI^l, we construct 
a genuine "gauge" invariant operator, which preserves the "gauge" invariance in the local observable universe. For the construc- 
tion, we note that the scalar curvature "i?, which transforms as a scalar quantity for spatial coordinate transformations, become 
genuinely "gauge" invariant, if we evaluate it in the geodesic normal coordinates span on each time slice. The geodesic normal 
coordinates are introduced by solving the spatial three-dimensional geodesic equation; 
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where ^T^^fc is the Christoffel symbol with respect to the three dimensional spatial metric on a constant time hypersurface and 
A is the affine parameter. We consider the three-dimensional geodesies whose affine parameter ranges from A = to 1 with the 
initial "velocity" given by 



dX 



^ g-C(A=0) 

A=0 



,-57(A=0)/2 



. (2.24) 



A point X' in the geodesic normal coordinates is identified with the end point of the geodesic, x'^{X, A = 1) in the original 
coordinates. Using the geodesic normal coordinates X\ we perturbatively expand as cc' = X' + Sx'^{X). Then, we can 
construct a genuinely "gauge" invariant variable as 

mix) 'R{t, x\X)) = X' + Sx'{X))) , (2.25) 

where t denotes the cosmological time and we introduced the abbreviated notation X := {t, X}. As long as the deviation from 
the FRW universe is kept perturbatively small, we can foliate the universe by the geodesic normal coordinates. 

Next, we quantize the system to calculate quantum correlation functions which become observable. A straightforward and fre- 
quently used way to preserve the gauge invariance is to eliminate gauge degrees of freedom by fixing the coordinates completely. 
In the local observable universe, the complete gauge fixing requires us to fix the boundary conditions in solving the constraint 
equations. However, to quantize the locally restricted system, we need to abandon several properties which are available in the 
quantization of the whole universe. One is that the quantization restricted to a local system cannot be compatible with the global 
translation symmetry at least in a manifest way. Another but related aspect is that basis functions for a mode decomposition 
become rather comphcated (even if it exists) than the Fourier modes. To preserve the global translation symmetry manifestly, 
we take another way of quantization than the quantization in the completely fixed gauge. 

In this paper, we perform quantization in the whole universe with the infinite spatial volume. The global translation symmetry 
in the spatial directions is then manifestly guaranteed in the sense that a shift of the spatial coordinates a; to a; + a just changes 
the overall phase factor in the Fourier mode by e'*'". Based on this idea, we consider and calculate observable quantities. Since 
the gauge invariant variable does not include the conjugate momentum of C,, we can consider products of at an equal 
time without the problem of operator ordering. The n-product of ^K, i.e., ^R{xi) ■ ■ ■ ^i?(a;„), preserves the "gauge" invariance 
in the local universe. To calculate the n-point functions of ^R, we need to specify the quantum state as well. One may think 
that the quantum state should be also invariant under the residual "gauge" transformations. However, we cannot directly discuss 
this invariance as a condition for allowed quantum states in this approach, because the residual "gauge" degrees of freedom are 
absent when we quantize fields in the whole universe. 

Here we note that even though the operator ^R is not affected by the residual "gauge" degrees of freedom, this does not imply 
that the n-point functions of f/? are uncorrected to the fields in the causally disconnected region. To explain this aspect more 
clearly, let us consider the n-point function of ^i? at i = tf whose vertexes are located within the observable region /. For a 
later use, we refer to the spacetime region which is causally connected to spacetime points in 0/ as the observable region 0. 
After expanding of the operator ^R in terms of (j, the interaction picture field of (, even if interaction vertexes which affect 
^R{tf, x) are confined within the observable region 0, the 7i-point functions of ^i? have correlation to the outside of through 
the Wightman function of C,i, which can be expressed in the Fourier space as 

G+(x-i, X2) - {Ci{t, x,)Ci{t, X2)) = I ■^e^''<"^^\vk{t)\^ , (2.26) 

where Vk is the mode function for Q. Since all the vertexes are confined within 0, the spatial distance \xi — X2\ is bounded 
from above. However, the IR modes with k < |a;i — a;2|~^are not suppressed and let the Wightman function diverge for 
scale-invariant or red-tilted spectrums. 

This long-range correlation becomes the origin of the IR divergence of loop contributions. Since the observable quantities 
should take a finite value, we request the IR regularity of observable fluctuations, which is achieved only when the quantum 
state is selected so that the long-range correlation is properly isolated from the observable quantities. In the following section, 
we will show that requesting the absence of the IR divergence in fact constrains the quantum state of the inflationary universe. 



III. IR REGULARITY CONDITION AND THE RESIDUAL "GAUGE" INVARIANCE 



A simple way to address the evolution of a non-linear system is to solve the Heisenberg equation perturbatively by assuming 
that the interaction is turned on at the initial time. In this section, taking this setup, we calculate the correlation function of ^R 
up to one-loop order and investigate the IR behaviour. 
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A. Specifying the iteration sclieme 



For notational convenience, we introduce the horizon flow functions as 



Pi. 
P 



1 d 

c 

£n-l dp 



(3.1) 



with n > 1, but without assuming these functions are small we leave the background inflation model unconstrained. The 
assumption that the interaction is turned on at the initial time requests 



CiU, X) = Cl{U, X) , TT{ti, X) = TTl{t„ X) 

where tt/ is the conjugate momentum defined from the non-interacting action: 



So = A/pi J dtd^x e'Pp'^ei 



(3.2) 



(3.3) 



as TT/ := 2M^^e^Pei(^i. Employing the initial conditions i3.2i . we can relate the Heisenberg picture field to the interaction 
picture field C/ as 



at,x) = Uj{t, U)Ci{t,x)Uiit, U). 



where Ui is the unitary operator given by 



Ui{ti, is) =Texp 



dtHi{t) 



(3.4) 



(3.5) 



Since the Heisenberg fields are related to the interaction picture fields by the unitary operator, the canonical commutation 
relations for the canonical variables ( and tt guarantee those for (j and ttj as 



[(lit, x), ^/(t, y)] = tS^^^x - y) , [C/(t, x), y)] = [7r/(t, x), mit, y)] = ■ 



(3.6) 



We note that the ri-point functions for ( given in Eq. ( 13.4b agree with those given by the solution of ( written in terms of the 
retarded Green function. We express the equation of motion for ^ schematically as 



where L denotes the differential operator 



(3.7) 



(3.8) 



and the left-hand side of Eq. ( 13. 7b is the same equation of motion as is derived from the non-interacting part of the action (13. 3b . 
By using the retarded Green function: 



Gr{x, x') = i9{t-t')[Ci{x), C/(.t')] , 



which satisfies 



(3.9) 



(3.10) 



with (5'^'*' {x — x') := S{t — t')5'-^'> {x — x'), the solution of the equation of motion with the initial condition ( 13.2b is given by 

C{x) = Ci{x) + L-^^2>^i^[x) (3.11) 

where the non-linear term is given by 

£^iSNL(i, x) := -2A/p2 ^ dt' j A^x'ei{t')e'P^''^p{t'fGR{x, x')§nl(x') . (3.12) 
Using Eq. ( 13.111 ). we can solve the equation of motion iteratively. 
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We expand the interaction picture field Q, which satisfies 

LCiix)=0, (3.13) 



as follows 



CKa;) - / 77^^ {a^Vke^'^-^ + h.c.) . (3.14) 



(27r)3/2 



The mode function Vk satisfies LkVk = where 



The mode function is normalized as 



g-2p 

Lk:=dl + {3-ei+e2)dp + ^k\ (3.15) 

{vke'''-^, Vpe'P-^) = {2TTf6^^\k-p), (3.16) 
where the Klein-Gordon inner product is defined by 

J d'xe^Pe,{f,dtf; - {dtfi)f;} • (3.17) 

With this normalization, we obtain the commutation relations for the creation and annihilation operators as 

[flfc, a^] =(5(3)(fc-p), [flfc, ap]=0. (3.18) 
Inserting Eq. (13.14b into Eq. ( I3.9l l, we can rewrite the retarded Green function as 

where Rk{t, t') is given by 

Rk{t,t') vk{t)vl{t') - vl{t)vu{t') . (3.20) 

We calculate the n-point function of ^R, setting the initial state to the vacuum defined by afc|0) = 0. 

Next, we explicitly calculate non-linear corrections. Since we are interested in the IR divergence, we pick up only the terms 
which can contribute to the IR divergences, employing the iteration scheme of For the convenience, we here introduce the 

symbol "w" as in 101 to denote an equality under the neglect of the terms unrelated to the IR divergences. The IR divergences 
mean the appearance of the factor {Q). For the scale invariant spectrum, this variance diverges logarithmically. Once a temporal 
or spatial differentiation acts on one of two C/s in the variance no longer diverges. In this sense we keep only terms which 
can yield {Q). When we write down the Heisenberg operator C in terms of the interaction picture field C/; at the one loop 
level, this is equivalent to keep only the terms without temporal and/or spatial differentiations and the terms containing only one 
interaction picture field with differentiation. 

In this paper, we refer to a term which does (not) contribute to the IR divergences simply as an IR (ir)relevant term. In the 
following, we will use 

TT? 

^'rCiDCi ~ Ci^rDCi , (3.21) 

where D is a derivative operator which suppresses the IR modes of ^/ such as dp and d"^ / pe^. Equation (|3.21t can be shown as 

• -1 , 



follows. The Fourier transformation of Lj^ CiDCi is proportional to 

d^p f dt'e,{t')e''>^'">p\t')Rk{t, t')Cip{t') {DCj)^_^ it') 



where C/fc and {DC,i)k denote the Fourier modes of C/ and DC,i. Since {DC,i)^,_p (f) — {DC,i)i^ (f) is suppressed and dp 
becomes time independent in the limit p ^ 0, the IR relevant piece of the integrand of the momentum integral can be replaced 
by 

Cip f dt'ei{t')e^P^''^p^t')Rk{t, t') {DCi\ (f) . 
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To obtain Eq. ( I3.21l i. we will also use 



£-V(x)«0 for/(x-)«0. 



Keeping only the IR relevant terms, the ADM action is simply given by 



{dpi? 



-2(p+C) 



where we used the lapse function N and the shift vector Ni, given by solving the constraint equations as follows. 



1 TT) p-2p 

p 



(3.22) 
(3.23) 

(3.24) 
(3.25) 



The action ( 13.23b can be easily derived from the action for the non-interacting theory. Extending the action ( 13.3b to a nun-linear 
expression to recover the dilatation symmetry, we obtain 



e-2(p+C) 



(3.26) 



where the abbreviated terms do not appear in the non-interacting action, which means that they should have more than two fields 
with differentiation. Such terms are IR irrelevant up to the one-loop order Thus, we obtain Eq. (13.23b . 
Taking the variation of the action with respect to C,, we obtain the evolution equation as 

By expanding C C ~ Ci + C2 + d ' ' ' ^ the equation of motion is recast into 

£C/=0, 

IR 



£C2 « -2C/AC/ , 



IR 



£C3 « -2 (C2AC7 + C7AC2 - CfAC/) 



where we defined 



(3.27) 
(3.28) 
(3.29) 



Solving Eqs. ( ll28T l and ( [l29l ), we obtain 



V. —d, , A := V,V, . 



C2 « -2C/£ .^AC/ 



TR 

Ca « C/^h' (4A£^iA + 2A) Ci 
where we noted the properties of given in Eqs. (13. 21b and (13. 22b . 



(3.30) 

(3.31) 
(3.32) 



B. Calculating the "gauge" invariant operator 

Solving the three dimensional geodesic equations, we obtain the relation between the global coordinates and the geodesic 
normal coordinates X*as 



Using the geodesic normal coordinates, the gauge-invariant curvature '-'R is expressed as 



IR 



-C(t.e-«X) 



X). 



(3.33) 



(3.34) 
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Here, we introduce the spatial average of ( at the local observable region whose spatial scale is Lt in the geodesic normal 
coordinates as 

Jd'XWL,{X)at,e-^X) 

^^'^ ■= — Wxmjx) — ' ^^-^^^ 

where Wl^ (X) is the window function which specifies our observable region. We approximate this averaging scale at each time 
t by the horizon scale, i.e., Lt ~ {eP^^'f p{t)}^^ . Using ^C(^)' we decompose the spatial coordinates as 

9i?(t, X) S m{t, e-{''"(*-'="'^)-'f(*)>e-''fWX) . (3.36) 
Using the window function Wl^ (X), we define a local average of an operator 0{t, X) as 

Id^xw,Ax)0{t,x) 

We note that the factor e"^^'*''^ ''X)-%{t)}^ which can be expanded in perturbation as 

^-{C{t,e-<X)-%{t)} 

= 1 - {(C/ - C/) + (C2 - C2) - iCiX ■ dxCi - CiX ■ dxCi} + ^(C/ - C/)' + • • • , (3.38) 

yields only IR regular contributions up to the one-loop order such as {{(j — C/)C/) ™d {CiX ■ dxCi)- The latter one also 
becomes finite because of the derivative suppression and the finiteness of | X | ' . For the reason mentioned above, neglecting the 
factor e"''^^^*''^ ^^)-%it)}^ vve rewrite the gauge invariant operator as 

m{t, X) « 'R{t, e-'f(*)X) . (3.39) 



Since the spatial curvature '^R is given by 



'R{x) « -Ae-'Pe-'^d'C , 



using the curvature perturbation in the geodesic normal coordinates: 

X{t, X) C(t, e-'^^^'^X) , (3.40) 
we can describe the gauge-invariant spatial curvature ^R as 

TO 

« -Ae-^Pd''x%{i. X) (3.41) 

at least up to the third order in perturbation. At the second equality, we again used the fact that the exponential factor e"^^'' ('^'^ 
do not give IR relevant contributions. As is pointed out by Tsamis and Woodard in Ref. iIstIi and also in Ref. | f38l |. using the 
geodesic normal coordinates can introduce an additional origin of UV divergence. We suppose that replacing C, in the geodesic 
normal coordinates with ^Q, which is smoothed by the window function, can moderate the singular behaviour. 
Inserting the solution of Q given in Eqs. (13.311 1 and ( 13.32b into Eq. ( 13.401 ). we obtain 

9C2«-C/(2£fl'A + X.ax)C/, (3.42) 

TR 1 — 9 — 

-Cj{2L-j,^/^ + Xdx) C/ + Cl(l-'C]^i£)X.ax£]^'AC/. (3.43) 
In deriving the expression of we used 

[£, = -2A. (3.44) 

We can verify that the term with (1 — L^L) in Eq. ( 13.431 1 vanishes. Since this term is eliminated by operating £, it always 
vanishes if its value and its first time derivative are both zero at the initial time, which is automatically satisfied by the definition 
of the retarded integral. 



One may think that interaction vertexes located outside the past lightcone can appear from the inverse Laplacian , which are included in the lapse function 
and shift vectors. In that case, |.X"| becomes unbounded. However, tuning the boundary condition of the constraint equations, which are the elliptic type 
equations, we can manifestly shut off the influences from outside the observable region 0. Since this change can be though of as the change of the residual 
coordinate degrees of freedom, it does not influence on the "gauge" invariant operator ^R. 
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C. Conditions for the absence of IR divergence 

We next discuss the conditions that the IR divergence does not arise in the expectation values of the gauge invariant variable 
^R. Hereafter, since we calculate the fluctuation only in the geodesic normal coordinates, we use x to express the geodesic 
normal coordinates. Then, using the above expression, two point function of up to the one loop order is obtained as 

TR — / 

{mixirRix2)) « 8e-4''(C|)9(\)a('2)(2(2£^iA(i) +a;r9,,)C/(.^i) (2£^1A(2) + a;2-5,.)C/(^2) 

+ (2£^;A(i) + xrd^yO{xi)Ci{x2) + Ciixi) (2£;^/A(2) + x^d^^fdix^)), (3.45) 

where we used an abbreviated notation d'^^-^ := ff^^ and Ajaj := {s'^p)~^d^^ for a = 1,2. As we mentioned at the end of the 
preceding section, the correlation function of ^R can contain IR divergences. One may think that the IR regularity is guaranteed 
if {2L^^A + xdx)0 = is imposed. However, this condition is in conflict with the use of the ordinary retarded integral in the 
iteration process. In fact, one can calculate (2£^^ A + x-dx)Qi as 

{2L-^'A + x-dx)Ci^ I ^^{e*-a.(^-2£^i,-^ + z,a..fc).,+h.c.| , (3.46) 

where -C^^^j. is the Fourier mode of The requirement that this expression should identically vanish leads to 

-2£fl'fc;^+*^-^)«fc = 0- (3.47) 

The first term is independent of x, while the second one manifestly depends on x. This shows that the condition {2L~^^A + x- 
dx)Ci = is incompatible with the use of the retarded integral. 

In place of this naive condition, an alternative possibility one can think of is to impose 



j^-^ {akDe'^-^Vk + h.c.) 



(3.48) 



(27r)3 
with 

D := /c-3/2g-»0(fc)fc . 5^fc3/2g#(fe) ^ (3 

where (f>{k) is an arbitrary real function. This condition continues to hold once it and its time derivative are satisfied on a certain 
time slice because both sides of this equation vanish under an operation of the second order differential operator L . The condition 
can be rewritten as a condition on mode functions 

-^^B'k-^^'k-Dvk, (3.50) 

which is not contradictory as in the case of Eq. (13.47b . With these conditions, the expectation value in Eq. ( 13.45b can be 
summarized in total derivative form as 

{m{x,rR{x,))'S8e^'''{-Cj) J ^^^i|^a,^,,{fc>.|V--(---)} , 

where we used k^^^D ~ e~"^'-'^''fc • (9fcfc'^/^e"*('^^ and J dfifc denotes the integration over the angular directions of k. Since this 
integral of total derivative vanishes, the IR divergence disappears. 

Although requesting the condition ( 13.50b can make the IR divergence vanish, a bit more careful thought rules out this 
possibility. Since the left hand side of ( 13.501 ) and its time derivative vanish at the initial time, the condition ( 13. 50b requests 
Dvk{ti) = Dvkiti) = 0. Operating the differentiation k ■ dk on the normalization condition of the mode functions, 

-2ie''Peiivkvl~Vkvl) = l (3.51) 

leads to a contradiction. The right hand side trivially vanishes after the operation of k ■ dk, while the left hand side gives 3. 

Even though the IR regularity condition ( 13.50b cannot be compatible with the initial conditions ( 13.2b . it is still instructive to 
give an alternative interpretation of the condition ( 13.501 ). We adopted the initial condition ( 13.21 ) for (/n, which identifies the 
Heisenberg fields with the corresponding interaction picture fields at the initial time and selected the vacuum state for the free 
field at the initial time. We denote a set of operations which specify the quantum state as mentioned above by an iteration scheme. 



11 



Then, when we take the same iteration scheme in the C,/tt system, it is not obvious whether the vacuum states in two respective 
systems are equivalent or not. Before closing this section, we show that the condition ( |3.50t is identical to the condition that 
these two vacua are equivalent. 

Since the transformation from C/tt to (^/tt is a canonical transformation, the corresponding correlation functions for the 
physically equivalent quantum state do not depend on the choice of variables, i.e., 

(C(a:i)C(.T2)) = (C(t, e'x\)~at, e'x\)) . (3.52) 

To employ the same iteration scheme, we request that both of Q and C, are solved by using L~j^^ by identifying the Heisenberg 
fields with the interaction picture fields at the initial time. We expand the interaction picture fields for Q and which we denote 
as Q and Q, respectively, in terms of the same mode function Vk as 

C/(^) = y (^^"fc^'^^'''" + l^-c- ' (3.53) 

C/(^) - j J^^'^kvue"' '' + h.c. . (3.54) 

Further, we select the vacuum states that are erased under the action of and in the respective systems. We denote these 
vacuum states specified by and 5^ as | ) and | ) , respectively. Now we show that the equivalence between the two point 
functions, i.e., 

(0|C(xi)C(rE2)|0) = (6|C(t, e^a;i)C(t, e^a;2)|0), (3.55) 
yields the condition ( I3.50l l. We expand C,{t, e*a;*) as 

C(i, e'x) = C{x) + sx ■ d^l{x) + ^s\x ■ d^fax) + 0(s3) 

IR. 

« C/ - C/2-Cfl,'AC/ + s(2£^iA + X ■ d^Xi + ■■■ , (3.56) 
taking into account that the interaction Hamiltonian is shifted by —s. Then, the right-hand side of Eq. (13.55b gives 

(0|C(t, e^a;i)C(t, 6^0^2)1 0) 

S (0 |{C/(a:i) - Ci{xi)2L-„^A^,^Ci{xi) + ■■■ }{Ci{x2) ~ C/(a:2)2£^iA(2)C/(.T2) + • • • }| 6) 

+ s(0 |(2£^iA(i) + x^ ■ d^Xi{x^)li{y) + • • • | 0) + ,s( |6(2;)(2£^i A + y ■ dy)li{y) + • • • | 6) 

+ O(s^) , (3.57) 

where "• • • " denotes higher order terms in perturbation. Since the terms in the second line of Eq. ( 13.571 ) agree with the left-hand 
side of Eq. ( 13.55b . the other terms on the right-hand side of Eq. ( 13.57b should vanish to satisfy Eq. ( 13.551 ). The remaining terms 
of order s at the leading order in in Eq. ( 13.571 ) is given by 

^ [«^e-''=-H2£H'A(i) .9.,Ke''^-^ +«fce''^-U2£fl'A(2)+a;2 •a.J«;!e-^'^-=] 



(27r)3 



(3.58) 



where we performed integration by parts. Now it is clear that requesting Eq. (13.521 ) gives the same condition as requesting the 
IR regularity ( 13.50b . except for the irrelevant fc-dependence of the phase of the mode functions. 



IV. SUMMARY AND DISCUSSIONS 



We have focused on observable quantities which are invariant under the residual coordinate degrees of freedom. These degrees 
of freedom are left unfixed in the ordinary discussion of cosmological perturbation. The key issue for the correlation functions 
of the "gauge" invariant operator to be IR regular is shutting off the long range correlation between observable quantities and 
the fluctuation outside the observable region. This is the matter of how we choose the initial quantum state. The leading 
effect of the long range correlation on local observables is the constant shift of the so-called curvature perturbation (= the trace 
part of spatial metric perturbation). The constant shift of the curvature perturbation can be absorbed by the scaling of spatial 
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coordinates. Assuming that the interactions are shut off before the initial time, we investigated the conditions that guarantee 
the equivalence between two systems mutually related by this scaling, and found that these conditions also guarantee the IR 
regularity of observable quantities. 

We also found that these conditions are not compatible with the prescription where the vacuum state is set to the one for the 
corresponding free theory at an initial time. In this setup, the initial time is very particular time, at which the Heisenberg picture 
fields agree with the interaction picture fields. Then, the horizon scale at the initial time becomes distinguishable from other 
scales. It is, therefore, natural that the presence of the particular initial time breaks the invariance under the scale transformation. 
One possibility to avoid breaking the symmetry is sending the initial time to the infinite past. In fact, in the standard cosmological 
perturbation by default we adopt the Euclidean vacuum which is selected by sending the initial time to infinity with the ie 
prescription. In the forthcoming paper, we will discuss when we choose the Euclidean vacuum, whether the IR contributions of 
the n-point functions for the coordinate independent operator ^i? becomes regular or not. 
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